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Abstract 

Neurons have the capability of transforming information from a digital signal at the dendrites of the presynaptic termi¬ 
nal to an analogous wave at the synaptic cleft and back to a digital pulse when they achieve the required voltage for the 
generation of an action potential at the postsynaptic neuron. The main question of this research is what processes are 
generating the oscillatory wave signal at the synaptic cleft and what is the best model for this phenomenon. Here, it is 
proposed a model of the synapse as an oscillatory system capable of synchronization taking into account conservation 
of information and consequently of frequency at the interior of the synaptic cleft. Trains of action potentials certainly 
encode and transmit information along the nervous system but most of the time neurons are not transmitting action 
potentials, 99 percent of their time neurons are in the sub threshold regime were only small signals without the energy 
to emanate an action potential are carrying the majority of information. The proposed model for a synapse, smooths 
the train of action potential and keeps its frequency. Synapses are presented here as a system composed of an input 
wave that is transformed through interferometry. The collective synaptic interference pattern of waves will reflect the 
points of maximum amplitude for the density wave synaptic function were the location of the "particle" in our case 
action potential, has its highest probability. 
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1. Introduction 

1.1. Synapses 

The structure of a neuron can be divided in three 
parts, receptor (dendritic tree), effector (axon) and 
nucleus where the receptor and effector converge. The 
dendritic membrane forms synapses with the axon’s 
tips of other neurons, with the special characteristic 
that there is no cytoplasmic bridge between them. The 
question at this point is: how information is conserved 
in this discontinuity? It is wildly known that dendrites 
receive input from hundreds of axon tips of other 
neurons, combine the input, which is delivered to the 
axon [Pur] ((2008}. However, how dendrites combine and 
deliver this information is not know with certainty, then 
is function of the axon to transmit the output of the 
dendrites to other parts of the nervous system. 

The nature of information at the synaptic cleft 
and along the axons is very different in amplitude 
and time scale, however it is finely tuned in order to 
keep the flow and coherence of information. Dendrites 
from the postsynaptic neuron, receive an almost con¬ 
tinuous oscillatory wave of ionic current input from 


the connecting synaptic cleft that joins it with the 
presynaptic terminal and convert it to a discrete voltage 
pulse that travels along the axon to be converted again 
in an oscillatory wave |Bullock| ( |1993) >. The wave to 
pulse basis transform takes place at the interior of each 
synapse, the mechanism is here proposed as follows: a 
pulse arrives on the axon terminal, this energy allows 
the entrance of positive calcium ions, which are going 
to move vesicles that carry neurotransmitters. This 
vesicles release their content outside the neuron where 
the oscillatory periodic wave takes the form of a field 
of ionic current. The mathematical details of this trans¬ 
formation are presented in the following sections. The 
output of the dendrites is the sum of waves resulting 
from all synaptic cleft activity, which are then delivered 
to the initial segment of the axon |Hodgkin and Huxley | 
( | 1990) 1, it is also presented a mathematical approach to 
this waves convergence. 

There are two types of activity encoded at synapses, 
excitatory and inhibitory, depending on the type of neu¬ 
rotransmitter that is released in the synaptic cleft, this 
chemical wave if excite can let to synchronous behavior 
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and generate a positive voltage or in case of inhibition 
can let to asynchronous behavior and generate a nega¬ 
tive voltage in the postsynaptic neuron. It is proposed 
that at the synaptic cleft the signals are smoothed and 
convolved in order to create the oscillatory wave, that 
has the same frequency as the originating train but have 
the important property that can be synchronized with 
other synapses for amplification of information. This 
information is then transmitted in space and delayed 
at the initial segment of the axon |Pur| ( j2QQ8| ). Once 
the wave has acquired enough amount of amplification 
after synchronization, the axon responds to this wave 
input by generating a pulse train, where each pulse has 
the same amplitude. But as a train of pulses, they keep 
the information flow frequency. 


of the signal. 

The wave to pulse conversion process, depends on 
the steady state level of the total dendritic synchroniza¬ 
tion, and its coincidence with the refractory period of 
the pulse (about 1 msec), during which no amount of 
stimuli can induce another pulse. Subsequently, there 
is a period lasting several milliseconds, in which an 
additional stimuli can induce a second pulse, only if 
the second stimuli is larger than the first. Given all this 
constraints, the wave to pulse transformation is non¬ 
linear and variable in time, if a steady above-threshold 
current is passed across the membrane, the train pulse 
has a high frequency, but this frequency depend on 
many other parameters that make it variable in time. 


1.2. Information flow 

The effect of all synapses working together, can 
excite or inhibit the production of a pulse in the receiver 
neuron. As mentioned before, this two effects are 
called excitatory synapses and inhibitory synapses, 
depending on the averaging effect of all the activity, 
the general output can be excitatory or inhibitory but 
not both. For the purposes of this paper, this general 
inhibitory or excitatory state, is the result of destructive 
or constructive wave pattern synchronization. 


Until now, it is considered that if two or more 
small inputs are given simultaneously, the responses are 
simply added, and the system is said to be linear. If two 
otherwise identical inputs are separated in time, and if 
the responses are identical except for time of onset, the 
system is time invariant Eccles (1964|>. Nonetheless, 


the mechanisms of synchronization, its variables, and 
limiting factors are unknown. On the other hand, it 
must be taken into account that the coupling medium 
is going to be crucial in the synchronization process, 
in the case of the synapse is mainly water, given that 
depending on the nonlinear nature of the synapse, 
it is possible to have in phase summation between 
presynaptic and postsynaptic neurons even when they 
are in anti phase, or it is possible to have anti phase 
behavior between the two cells even when they are in 
phase. 


Another important effect that must be taken into 
account, is the secondary effect that the neurotrans¬ 
mitter release of a presynaptic terminal can have in 
neighbouring postsynaptic dendrites. There is a high 
possibility of having this cascade amplification given 
the analogous dispersive nature of the synapse signal 
and more importantly, given the required amplification 


1.3. Cable Equation and chemical synapses 

The cable equation is a one-dimensional diffusion 
equation, and it has been widely used to model the 
dynamics of diffusion at chemical synapses. When 
a pulse discharges at the axon terminal, some of the 
vesicles in the pre-synaptic neuron discharge their 
contents into the synaptic cleft, this discharge has been 
understand as homogeneous. In the following sections 
I explain the oscillatory nature of this discharge, and in 
the characteristics that can be observed in the coupling 
between oscillators. 

There is no biological proof of this proposed the¬ 
ory because the chemicals, ions and neurotransmitters 
directly in the synaptic cleft have not been measured. 
However, the concentration of this elements have been 
measured at the post-synaptic terminal, the concen¬ 
tration rises rapidly to a maximum value and then 
decays slowl}{Pur] ( |20Q8) >, components of the chemical 
information can be detected many milliseconds after a 
maximum of activity. 

1.4. Synchronization of synapses 

Connectivity among neurons is required, and is 
present in massive amounts in order to generate 
synchronous coupling and amplification, with enough 
amplitude to transmit the analog oscillatory information 
of chemical synapses |Hebb| ( fl949| l. This activity should 
be nonlinear and self-sustaining, as was observed in 
biological oscillators. 

This nonlinear synchronization behavior was first 
studied in a bigger scale than in this thesis, and 
only taking into account the diffusion coupling by 
Katchalsky |Katchalsky| (fT976a]l. He pointed out that a 
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large number of nonlinear elements diffusely coupled, 
will give their properties of energy flow to the whole 
system. In this way, the system stabilises in a zero 
equilibrium state for E=0. Subsequently, the system 
is driven progressively away from equilibrium, until a 
critical level is reached where a phase transition occurs 
for E + 0. 


2. Steps of the theory 

Proposal and modelling of synapses as non linear os¬ 
cillators capable of synchronization. 

• Research on the possibility and constraints for os¬ 
cillation in the synapse. 

• Oscillator model of a synapse 

• Propose coupling between synapse 

• Synchronization constraints along synapses 

• Relation between spectrum of the total synaptic 
signal and the frequency spectrum of the action po¬ 
tential. 


2.1. Synchronization in the neri’ous System 

Biological oscillators like neurons and heart cells are 
usually modelled as nonlinear systems in a dimension 
2 space. It is possible for this models to present limit 
cycles and synchronization. Nonlinear electronics as 
the Van Der Pol Oscillator and the Josephson function 


are good models for neuronal oscillators (Strogatz 

12001 ) . 


A neuron has cycles of activity that oscillate be¬ 
tween states of above threshold and sub threshold 
activity. Once the cell has the required energy, at a 
point in its cycle it releases an electrochemical signal 
that enables the communication channel to connect 
with other neurons in a synchronous fashion. Using the 
phase space of the neuron, let the description of how 
close the neuron is to reach the above threshold firing 
state. This cycle has a characteristic period and an 
amplitude. This means that the oscillations will remain 
near a constant frequency, even when the neuron is 
perturbed (Matthews et al.|[l991[l. 


rise to the field of nonlinear dynamics given that he ob¬ 
served that all initial conditions of the triode generator 
converged to the same final periodic orbit. In trying to 
model this behavior, they found equations of the phe¬ 
nomenon that couldn’t be solved as we are used to solve 
linear equations, nonlinear Van Der Pol system of equa¬ 
tions. 

X + /aX(X 2 - 1) + X = 0 (1) 

From the first and last term we can have an idea of the 
system that we are analyzing, and how we are going 
to model it because this two terms reflect the nature of 
an harmonic oscillator that can be achieved by the be¬ 
haviour of an inductance and a capacitance. However, 
the term in the middle reflects a different oscillatory be¬ 
havior, this is a damping nonlinear term that is origi¬ 
nated by the functioning of vacuum tubes. We can see 
that the damping in the system can be positive or nega¬ 
tive depending on the value of X. If X > 1, it is a normal 
damping that tends to make the amplitude of the oscil¬ 
lations decay. But, if A < 1 it is a negative damping 
which means pumping, this tends to amplify the ampli¬ 
tude of the oscillation. This pumping nonlinear effect is 
desirable for modelling the synapse given that neurons 
have a perpetual non zero activity that can be modelled 
through this term. 


3. Description of the biological phenomena 

The main interest of this research is to answer the 
wave to pulse generation problem in the chemical 
synapses of our nervous system. The current state 
of the art in this respect is pretty scarce and unclear, 
regarding conservation of information and frequency at 
the interior of the synaptic cleft. The curiosity to solve 
this problem was mainly raised by the fact that the 
trains of action potentials certainly encode and transmit 
information along the nervous system but most of the 
time neurons are not transmitting action potentials, 
99 percent of their time they are in the sub threshold 
domain were only small signals without the energy to 
emanate an action potential are the ones that carry the 
majority of information, the one that let us perceive 
the world in one way, the same synchronised way 
that let us have a language, memory and in general, 
activities that do not require the fast response inter 
neuron communication as in electric synapses. 


2.2. First theoretical study of Synchronization. 

Balthazar van Der Pol studied Synchronization theo¬ 
retically for the first time. He worked synchronising tri¬ 
ode generators from vacuum tubes. This research gave 


The model that is proposed in this thesis for a 
synapse, is an oscillatory diffusion of ionic current, 
that smooths the train of action potential and keeps 
its frequency. Synapses can be seen as nonlinear. 
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self-sustained oscillators with stable limit cycles, the 
linear oscillator, which can cycle at any amplitude, is 
not in agreement with real biological oscillators that 
have a regulated amplitude. The network of synapses 
should synchronise to generate the Gaussian properties 
of synchronous oscillators, this is achieved after the 
realisation of a Fourier series that transform from the 
analogous input of information that is given in each 
synapse oscillation to the almost discrete response 
along the axon in the form of an action potential. 
Synapses with the axonsAZs tips of other neurons, with 
the special characteristic that there is no cytoplasmic 
bridge between them. The question at this point is: how 
information is conserved in this discontinuity?. 

The nature of information in dendrites and axons 
is completely opposite, yet the same. Dendrites 
receive an almost discrete pulse input and convert 
it to an analogous continuously oscillatory wave of 
ionic currcn (Bullock| ( |1993| ). The wave to pulse basis 
transform takes place at the interior of each synapse, 
the mechanism is as follows: a pulse arrives on the axon 
terminal, this energy allows the entrance of positive 
calcium ions, which are going to move vesicles that 
carry neurotransmitters. This vesicles release their con¬ 
tent outside the neuron where the oscillatory periodic 
wave takes the form of a field of ionic current. The 
output of the dendrites is the sum of waves resulting 
from all pulse inputs, which is delivered to the initial 
segment of the axon |Bullock| ( fl993] > . Locally, at the 
synapses, the signals are smoothed and convolved in 
order to create the oscillatory wave, that has the same 
frequency as the train but can be synchronized with 
other synapses for amplification of information. This 
informa- tion is then transmitted in space and delayed 
at the initial segment of the axon [Pur 2008]. Once the 
wave has acquired the enough amount of amplification 
after synchronization, the axon responds to this wave 
input by generating a pulse train, where each pulse has 
the same amplitude. But as a train of pulses, they keep 
the frequency information flow. 

The purpose of this thesis is the explanation and 
demonstration, that synapses synchronization is nec¬ 
essary to generate a train of pulses along the axon. 
Addition- ally, this synchronization is achieved when 
the spectrum of all the signal from the synapses, 
matches the frequency spectrum of the action potential. 
Presynapse as oscillator postsynapse as oscillator 
and synapse as medium that couple them to create 
synthronization through the transmission of motion. 


A neuron oscillates through oscillations in its elec¬ 
trical field, and at a point in this cycle it releases an 
electrochemical signal to connected neurons. The 
phase space of the neuron describes how close the 
neuron is to firing this signal, and this cycle has both a 
characteristic period and an amplitude. This means that 
the oscillations will remain near a constant frequency, 
even when the neuron is perturbed. 


4. Theretical proposal 


Here I analyse the case in which a train of action 
potentials in the presynaptic neuron does not cause the 
postsynaptic one to fire, which is about 99% of the neu¬ 
ron’s life Katchalskyl ( 1976b} and ?. When a train of 
action potentials arrives at the synapse in the presynap¬ 
tic neuron, it is converted into a chemical signal and this 
chemical signal in turn is converted into a small electri¬ 
cal signal that preserves the original frequency from the 
action potentials, see figure ??. 



Figure 1: A train of action potentials in the presynaptic neuron is 
smoothed out by being converted into an oscillatory signal of neuro¬ 
transmitter concentration in the synaptic cleft. This neurotransmitter 
concentration is then turned into a low amplitude oscillatory electrical 
signal in the postsynaptic neuron. The original frequency of the train 
of action potentials is preserved. 

4.1. The synapse as a system 

When analyzing this process in more detail we 
observe that when the action potential reaches the tip 
of the axon of the presynaptic neuron, it triggers the 
movement of vesicles that release neurotransmitters 
that then travel through the synaptic cleft via a diffusive 
process towards the postsynaptic neuron. Afterwords 
an electrical signal is converted into a diffusive process 
modulated by the frequency of the pulse. Once the 
neurotransmitters reach the postsynaptic neuron, this 
oscillatory chemical signal is converted back into 
an oscillatory electrical signal. This process can be 
understood in the sense of chemical oscillations ?. 

This process can be modelled interpreting the train of 
action potentials as the input signal, and what happens 


4 


























with the vesicles, the synaptic cleft and the receptors 
can be interpreted as a system that transforms the input 
signal into a smoothed signal with lower amplitude. 
This idea can be seen in figure [2] 



SYSTEM: 

Vesicles, 

neurotransmitters, 
synaptic cleft and 
receptors. 

v___y 


Figure 2: Vesicles, neurotransmitters and synaptic cleft as a system 
that transforms the train of spikes into a smoothed out version with 
the same primary frequency. 


The following assumptions will apply for the model 
in figure [2] 





Neurotransmiters concentration pulse 
spreads as it travels through the synaptic cleft 


Figure 3: A neurotransmitter concentration leaves the presynaptic 
neuron with a certain velocity and as it travels through the synaptic 
cleft, it spreads by diffusion its concentration. 


Therefore, similar to what has been described above, 
the postsynaptic neuron receives a smoothed concentra¬ 
tion wave, which is then converted into a small periodic 
electrical signal, proportional to the original concentra¬ 
tion. 


• The post-synaptic neuron will not emit an action 
potential, but will transmit small subthreshold sig¬ 
nals. 


• The synaptic cleft dimensions stay the same, plas¬ 
ticity or changes in the configuration of the synapse 
are not happening. 


If the above conditions are met the system in figure 
[2]behaves as a time invariant system, which means that 
its output can be written as the input convolved with the 
impulse response of the system. In short words we can 
apply Fourier analysis to it. Thus, we can write: 


Output 

= 

Input 

* 


Impulse Response of the System 


( 2 ) 


To obtain the impulse response of the system, we 
need to analyse what happens at the synaptic cleft in 
more detail. When vesicles release neurotransmitters 
into the synaptic cleft, they release the concentration 
with approximately the same velocity at which the vesi¬ 
cles fused with the cellular membrane. As this blob 
of concentration travels through the synaptic cleft, it 
spreads by diffusion as seen in figure [3] 


4.1.1. Vesicles and Diffusion at the synaptic cleft 
The concentration gradient moves through the synap¬ 
tic cleft, and as it moves it also diffuses, remembering 
that the medium is 80% water at the interior of the 
cell. Thus, the postsynaptic neuron receives a diffused 
concentration potential. Assuming that the concentra¬ 
tion pulse is on a fluid, moving at constant velocity 
for simplicity, the dominant effect on the spread of the 
signal is diffusion. 

Diffusion in the direction along the synaptic cleft can 
be quantified, in the frame of the moving concentration 
pulse, by the 1-D diffusion equation: 

dc(x,t) d 2 c(x,t) 

dt ~ dx 2 1 ’ 

This equation is solved applying the Fourier transform 
in space using: 


T{g (n) (x)} = (2 nisTTg(s) 
for the derivative. Thus equation [3]becomes: 

= _2(4) 
dt 

Where ffc^s, t), is the Fourier transform of the concen¬ 
tration. Equation [4] is a first order differential equation 
in time, therefore: 

Tc(s, f) = Tc(s, 0)<T 27rVz5f (5) 

To solve this equation we use the following Fourier 
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Applying the inverse Fourier Transform to equation [5] 
we get: 


T~ l Tc(s,t) = r~ l Tc(s, 0) * T-V 21 ^ 2 - 0 '} 

1 


c(x, t ) = c(x, 0) * 


V2 nDt 


e iD, 


( 6 ) 


Equation [6] prescribes how the concentration pulse 
spreads as it travels through the synaptic cleft. For now 
we assume that the synaptic cleft has a defined size and 
therefore the concentration pulse travels a predefined 
distance, which we will call d c . Additionally the vesi¬ 
cles release the neurotransmitters at a velocity v„ that 
in average is a constant. Therefore, the concentration 
pulse spreads over a time t that we can predict: 


d c 

t — — ,d c = Synaptic cleft width, 

v„ 

v r i = Velocity of neurotransmitters wave 


The equations described in this model are a simplifica¬ 
tion of the main idea, that concentration pulses travel the 
synaptic cleft at a more or less standard time, since most 
synaptic clefts exhibit the same dimensions. Thus we 
can write the total spread, that the concentration pulse 
acquires while traveling through the synaptic cleft as 
follows: 


c(x,d c /v„) = c(x, 0) * 


1 =£ 
- - g 2U.dc In ,) 

y/2nD(d c /v„) 


(7) 


Additionally, we can now rename the diffusion constant 
D, d, and v n by one single constant, d, that quantifies 
the spread of the signal as it travels through the synaptic 
cleft, as shown in the following equation: 


c(x,d c /v „) = c(x, 0) * 



( 8 ) 


We have to remember that this solution is derived in 
the coordinate system of the pulse where the concen¬ 
tration pulse is emitted by the presynaptic neuron and 
is absorbed by the postsynaptic neuron. Furthermore, 
the concentration of neurotransmitters will be a func¬ 
tion of the voltage in the neuron where, c(t ) = A pit) 
and with this in mind, we can rewrite the concentra¬ 


tion pulse in the neuron as a voltage function of time 
for both the pre and postsynaptic neuron: The following 
equations describe the following sequence of events 1- 
Voltage converted to concentration pulse. Pre-synaptic. 
2-Concentration converted to voltage. Post-synaptic. 3- 
Spatial convolution converted to time convolution at the 
synapse. 


c(x,0) = A m p m (t), (9) 

c(x, d c /v„) = A out p out {f), (10) 

-v 2 -Min ') 2 

e 2 « = e 2 f , (11) 

As mentioned before, I want to reiterate that the neuron 
converts voltage signals to concentration signals at the 
presynaptic neuron and concentration signals to voltage 
at the postsynaptic neuron. This happens in a similar 
way, as a video camera in our phone converts a visual 
signal into a digital one, and then the screen converts a 
digital signal into a visual signal. Based on this analogy, 
the spatial convolution in equation [8] can be converted 
into a time convolution: 


AquiP 'out(f) ~ A In Pi nit) * 




( 12 ) 


We have differentiated the voltage-to-concentration and 
concentration-to-voltage conversions at the pre and 
post-synaptic neuron respectively, because these two 
types of conversions are governed by different under¬ 
lying processes. Rearranging equation[l2] we have: 


Pout(t) = -— Pin(t) * 
A out 



(13) 


Equation [13] explains what happens to the action po¬ 
tential as it reaches the end of the presynaptic neuron, 
which is converted to a chemical signal and then back 
into an electrical signal. The scaling constant 1 /A out 
quantifies how much the amplitude of the signal is re¬ 
duced. Additionally, the Gaussian function present in 
equation[l3]is normalized, thus the convolution smooths 
out the signal. We now have a mathematical represen¬ 
tation of the system in figure [2] To understand better 
what happens in equation [2] we can appeal to figure [4] 
where the input signal of spikes or pulses is smoothed 
and becomes small. The system box in figure [4] depicts 
the behaviour of the synapse. 
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Figure 4: The synapse as a energy conversion and smoothing system. 


We now have a mathematical explanation for the sys¬ 
tem in figure[2]and the conceptual equation[2| where we 
have explained the vesicles, neurotransmitters, synaptic 
cleft, and receptors as a system quantified in our final 
result in equation [13] 

4.1.2. Predictive model for a train of action potentials 
crossing the synaptic cleft just after conversion 
into neurotransmitters. 

In the previous section we have just derived what 
happens at the synapse and now we apply the result of 
equation [13] to predict what will happen to a train of 
action potentials. 


Thus: 


Pout(f) = - — frit) * 
A-out 




06) 


2 S(t - kT) 




H /A* 




k=—co 


^/Ar 


-e 
i 
n 

(17) 

_,2 

(18) 


And using the fact that 6(t - a) * /(f) = f(t - a), we 
obtain: 



1 -<r -kT) 2 

e 2(IA l 

out V 

^/Ai 


(19) 


First we need to write a train of action potentials in 
mathematical terms and for now I propose to model 
them as a train of Dirac’s delta functions, as this cap¬ 
tures their main characteristics, especially the fact that 
they arrive at a certain frequency at the end of the synap¬ 
tic neuron. This idea is illustrated in figure [3] 



► E S(t-kT) 

k— — oo 


Figure 5: A train of action potentials written as Dirac’s delta func¬ 
tions. The time interval between two action potentials is T. 


We can express the train of action potentials as a train 
of Dirac’s delta functions yet, in reality we will have a 
bounded train of action potentials. However, for now 
and without loss of generality, I will write it as an in¬ 
finite train of deltas as it is a well known function in 
Fourier analysis: 


Action potentials spaced T = ^ 6(t - kT) = frit) 

k——oo 

(14) 

With this formalism for the action potentials, frit), we 
can now see what will happen to the electric signal as 
it travels through the synapse. Using our result from 
equation [13] and equation[l4j we have: 


This is a surprising result, since a train of action 
potentials has become a train of smooth Gaussian 
functions, which perfectly resembles a wave. Now, 
what we quantitatively saw in figure [4] is expressed in 
equation [19] The synapse converts a train of action 
potentials into a smoothed oscillatory signal. From 
all of this follows that the behaviour of a subthreshold 
signal going into a postsynaptic neuron can easily be 
modeled by an oscillator that outputs a Gaussian like 
wave. Furthermore, this behaviour has been studied by 
expressing a train of action potentials as deltas ?. When 
we represent an action potential in more detail, the 
signal will be even closer to a wave and since a delta 
function is the sharpest signal that can exist, it can be 
turned into a smoother gaussian. Therefore, a train of 
action potentials, which is a continuous signal, will be 
more smoothed out by the synapse and the output will 
become closer to a sinusoidal wave. 

With this, I have demonstrated that for the purpose of 
studying the signals going into a single neuron, we can 
substitute the input synapse by an oscillator that outputs 
waves of a predefined frequency. 

5. Discussion 

Even though there is no biological proof regarding 
the synapse as an oscillator, there are many advantages 
for neural networks applications if we model the 
synapse as an oscillator. If we consider a detailed 
biophysical model of the synapse and we analyse the 



Pin(t) = frit) 


(15) 
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behavior of the different components, from calcium 
dynamics to vesicles release, it is possible to see the 
oscillatory activity. Those biophysical models make 
use of experimental data in order to solve the coupling 
of differential equations in a numerical way through 
the tuning of variables. In this way, we can be more 
confident about the assumption of oscillations along 
the presynaptic terminal, synaptic cleft and postsinaptic 
receptors. 

It is through wave interference patterns that this 
model generates the enough amount of fluctuations 
and coherence to represent the activity recorded in the 
postsynapse. 
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